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Abstract 


Introduction. Increasing the durability of gas turbine engine (GTE) blades is achieved through the use of special 
protective coatings on their surface. For the development of such coatings, the basic source information is the geometric 
profile of the blade section. To transfer a given blade cross-section profile to the appropriate CAD/CAM system or 
engineering analysis package, parametric modeling methods are used to automate this operation. However, the known 
approaches to creating a parametric model of a blade profile are not without a number of disadvantages, and a generally 
accepted method for creating it does not currently exist. The research was aimed at creating a technique for approximating 
the profile of gas turbine engine blades, convenient for use in the subsequent analysis of the operating conditions of 
special coatings on the surface of the blades. 

Materials and Methods. When constructing parametric models of the profile of gas turbine engine blades, a method based 
on the orthogonal Legendre polynomials was used. This made it possible to provide high accuracy of approximation and 
construction of a continuous mapping for the parameters of the blade profile approximation. A Python application was 
created for automated processing of source profiles. It provided the calculation of the coefficients of approximating 
polynomials for the contour lines of the blade, visualization of the calculation results, and creation of a dxf file based on 
the points of approximating functions to transfer it to the CAD system. Next, geometric models of blades were used to 
solve the problem of a stationary aerodynamic flow around a blade. The results of solving this problem were used to study 
the effect of the blade profile on its cooling in an aerodynamic flow. 

Results. As an example, three options of blade profiles belonging to different types of GTE were considered. It was shown 
that for all three studied profiles, the proposed technique provided obtaining parametric models that maintained high 
accuracy in constructing approximating lines, which was confirmed by the values of the determination coefficients close 
to unity. To illustrate the possibility of using the obtained models, examples of solving the gas dynamic problem with a 
potential flow around a blade in a stationary aerodynamic flow were given. The distributions of pressure and temperature 
on the surface of the blade were calculated using the finite element method. 

Discussion and Conclusion. The calculation results show that the proposed technique of approximating the profile of the 
GTE blade, based on the use of orthogonal polynomials, is a convenient tool to automate the creation of a geometric model 
of the blade and compare different types and profiles of blades, solving the corresponding gas dynamic problems. At the 
same time, for a given blade profile and GTE operating conditions, it is possible to obtain the distribution of temperatures 
and forces acting on the surface of the blade, which is required for predicting the durability of special coatings. 
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AHHOTalna 

Beedenue. Wospnuenue JONTOBeYHOCTH JOMaTOK ra30TypOHHHBIX ABuratTene (IT) qocTuraetca 3a cueT IpHMeHeHHA 
CHeUMaIbHBIX 3all[MTHbIX IOKPbITHH MX NoBepxHocTu. Jt pa3paOoTKH TaKHX MOKPbITHM OCHOBHOM HCXoOZHO 
MHopMallvet ABJIAeTCA TeOMeTPHYeCKHH Upodurs cevenua NoMaTKH. Ja Nepeqaun 3ayaHHoro Mpodusa cevenna 
louaTKH B cooTBeTcTByrollyIo CAD/CAM-cuctemy wim lakeT WHXKeHepHOrO aHasIM3a MIPHMCHAIOTCA MeTOIbI 
HapaMeTpH4ecKOrO MOJeJIMpOBaHHA, MO3BOJIAIOUIMe ABTOMATH3HpOBaTb JaHHy!O Orepaluio. OHako M3BeCTHBIC 
MOAXObI K CO31aHHIO MapaMeTpH4ecKON MOE NpOPUA WOMaTKH He JIMIMCHbI pia HeAOCTATKOB, H OOMIeNpHHATON 
MeTOAUKM ee CO3aHHa B HacTosuee BpemMa He cyujecTByeT. Llemb1o WaHHOM paOoTb! ABJIAeTCA CO3qaHHe MeCTOAMKH 
almpokcuMalun mpoduna sonatoK CTI, ynoOuol Wa Mcnob30BaHHA TPH MOCMeAyIOUeM aHasIM3e YCOBHM paOorTsl 
CIeMMasIbHBIX MOKPbITHH MOBepXHOCTH JIOMAaTOK. 

Mamepuanoi u memoooi. pu noctpoenuu MapaMeTpHyeckux MOJeseH MpodpusaA JOMAaTOK ras0TypOMHHBIX JBHraTesen 
aBTOPaMH MCHOIb30BaH MCTON, OCHOBaHHbIM Ha IPHMCHeCHHM OPTOTOHAJIbHBIX MOJMHOMOB JlexaHypa. ITO WO3BOIMIIO 
OOecIe4YHTb BbICOKYIO TOYHOCTS alllIpOKCHMallMH WM MOcTpoeHHe HelpepbiBHOrO oTOOpaxKeHHA JIA MapaMeTpoB 
allmpokcuMalun npopusa TomaTKH. J[1t aBTOMaTH3HpOBaHHOM OOpaooTKH HCXOAHbIX Ipoduse co3aHo IpH0xKeHHe 
Ha s3bIke Python, mosBosAIOlWee BEIMHCIATS KOOPPUUMEHTH! alMpOKCHMUPyIOWHX MOMMHOMOB JIA JIMHUM KOHTypa 
JIOMAaTKU, BH3YaIM3MpOBAaTb pe3yJIbTaTbI pacueta H CO31aBaTb 10 TOUKAM alllIpoKCHMUpylowMx dpyHKunl dxf-datin wa 
Tlepegqaun ero B CAD-cnctemy. JJaee reoMeTpHueckHe MOJCJIM JIOMATOK HCHONL3OBAaIIN JIA pewleHHaA 3aaun OOTeKAaHHA 
JOMWAaTKH CTalMOHapHbIM adposMHAaMHYeCKUM TOTOKOM. Pe3yJIbTaTbI PeWIeCHHA 3ZTOM 3aadH UCTIONb30BaHbI IIpv 
MCCIe{OBAHHH BIMAHHA IIpOpUsA JOMaTKH Ha ee OCTLIBAHHE B AaIPOMHAMMYECKOM HOTOKE. 

Pe3yibmamot uccnedosanua. B kadecTBe IpHMepa paccMOTpeHbI Tp BapHaHTa Mpopuset OMAaTOK, OTHOCALNUXCA K 
pasHpim TunamM ITJ]. Iloxa3aHo, 4ro fd BCeX TpeX H3Y4eHHBIX Mpodusei mpeqOKeHHAA MeTOAMKa MO3BOIAeT 
TOy4aTb MapaMeTpHyeckHe MOJeIH, OOeCHeYMBaIOMIMe BbICOKYIO TOUHOCTS MOCTPOeHHA alINIpOKCHMUpYIOWINX JIMA, 
4TO NOATBepxRyaeTcA OH3KHMH K eAMHUWe 3HaxyeHHAMM KOSdduuMeHTOB WeTepMuHauuH. JIA wisrOcTpayHn 
BO3MOXKHOCTH HCIHOJIb30BaHHA MOJYYCHHBIX MOJeNeH MpHBeAeHbI IpHMepb! peuleHHA 3aa4H ra30BOl AHHaMUKU TIpu 
TMOTCHIMaIbHOM OOTeKaHHH JIOMATKH B CTalMOHapHOM a3poyqMHamMMyeckom noTOKe. MeToOM KOHeYHBIX 3JIEMeCHTOB 
paccunTaHbl paciipeyeeHuaA JaBIeHHi M TeMiepatyp Ha MOBepXHOCTH JIOMAaTKH. 

O6cyotcoenue u 3aksu0o4enue. PesyibTaTbl BEIMMCICHHH NOKa3asIH, YTO TIpewiaraemMasd MeTOAMKa allmpoKCHMallMH Mpopwa 
nonatku IT], ocHopaHHat Ha HCMOb30BaHHH OPTOTOH@JIbBHbIX MHOFOYWICHOB, ABJIAeTCAH yOOHbIM HHCTpyMeHToM, 
TIO3BOJIAIOLIMM ABTOMATH3HpOBaTb CO31aHe reOMeTPHYECKOM MOJ{eJIH JIOMATKH H IpOBOAUTb cpaBHeHHe pa3JIM4YHbIX THIOB 
MW Mpoduset JoMaTOK, peliiat COOTBeETCTBYIOMIHe 3aayH Ta30Bol AHHaMUKU. IIpu 3TOM AWA 3ayjaHHOro UpopusA OMAaTOK 
MW ycnoBuli padots! IT] MoxkHo noslyuuTb paciipeyeseHue Temiepatyp UM ycuIMi, JevicTByIOWIMX Ha MOBeEpXHOCTH JIoMaTKH, 
4TO HeEOOXOAHMO JVI NPOTHOSHpOBaHHA JOJITOBCYHOCTH CIeLMaIbHBIX MOKPBITHH. 


Ksnrouesble CJI0Ba: jOnaTKH Tra30TYPpOHHHBIX Buratelen, pou ceqeHHaA, allupOKCHMalHaA, OpTOrOHaJIbHble 
MHOrouWJICHbI 


BaarogxapHocru. ABTODBI BbIpaxKaroT OaroapHOcTB peak U peweH3eHTaM 3a BHHMAaTeCJIBHOeC OTHOIMICHHE K CTAaTbe 
HW YKa3aHHble 3aMeuaHhA, YCTpaHeCHHe KOTOPBIX MO3BOJIHJIO MOBbICHTb Ce KadeCTBO. 


Aaa WaTnupoBanna. Conosses M.E., lynesa 1O.H., bangaes C.JI., bangaes JI.X. Anmpokcumana npopusa slomaToK 
Ta30TypOHHHBIX jBuraTelelti. Advanced Engineering Research  (Rostov-on-Don). 2024;24(1):78—-87. 
https://doi.org/10.23947/2687-1653-2024-24- 1-78-87 


Introduction. Turbine blades are the most highly loaded parts of a gas turbine engine (GTE), the durability of which 
mainly determines the engine overhaul time [1]. One of the ways to increase the durability of turbine blades is to apply 
special protective coatings on their surfaces [2]. The coating is a complex composite structure. To select materials for it 
and determine its optimal geometric characteristics, it is necessary to know the operating conditions of the part: the 
distribution of temperatures, pressures, and shear stresses on the surface of the blade. 
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The highest requirements are posed on the shape of the blades, as well as on the selection of materials, their 
manufacturing technique, and special coatings [3]. When designing special coatings, the basic initial information is the 
geometric profile of the surface of the part, since the parameters of the coating manufacturing process are set on its basis. 

Among the characteristic elements of the profile, the following can be distinguished: chord — distance between the 
most distant points with a horizontal arrangement of the blade, like an airplane wing; suction side — upper part of the 
profile; pressure side — lower part of the profile; leading and trailing edges [4]. 

The shape of the blade cross-section profile is selected both on the basis of experimental studies of pilot units and 
actual engines [5], as well as on the results of numerical experiments on models [6]. In the course of numerical simulation, 
we study the task of creating a geometric model of the blade based on a given cross-section profile and transferring it to 
the appropriate CAD/CAM system [7], or an engineering analysis package [8], in which the parameters of the physical 
properties of the materials of the parts are set, the boundary conditions of the problem are determined, and calculations 
are performed. 

Traditionally, when creating a blade drawing, its profile is described by a set of circular arcs, which is then transferred 
point-by-point to the module for developing a geometric model and generating a grid of finite elements. The disadvantage 
of this approach is that when building each new model, you have to manually create its profile in the graphical editor. In 
this regard, methods of parametric modeling of the blade cross-section profile are proposed to automate this operation. 
One of the most well-known models is the nine-parameter RATD (Rapid Axial Turbine Design) model described in [9]. 
At the same time, the practice of using RATD has revealed certain disadvantages inherent in this model, specifically, the 
inconvenience of using it when optimizing the geometry of the profile, as well as insufficient accuracy. On this basis, the 
authors [10] proposed a modified version of this model, which included a larger number of parameters. Methods of 
parametric modeling of blade cross-section profiles are also actively used in domestic practice. Thus, in [11, 12], 
algorithms for the automated construction of the blade profile using approximation by curves from a set of parabolas and 
Bezier curves of the second order were proposed. In [13], a method for designing a grid of GTE profiles based on these 
algorithms was described. The use of parametric models of blade profiles in solving optimization problems was 
considered by the authors in [14, 15]. 

At the same time, it should be noted that the existing methods of parameterization of the blade profile are not without 
a number of disadvantages. The main disadvantage of traditional blade profile parameterization schemes based on the use 
of second-order curves is that different functions have to be used for different parts of the workpiece to maintain the 
accuracy of the description. At the same time, these functions do not form an orthogonal system, and the values of the 
parameters determined by the approximation of the existing profile by the least squares method turn out to be correlated. 
This leads to the fact that it is not possible to construct a continuous mapping for a set of parameters when the system 
approximates the cross sections of the three-dimensional surface of the blade of double curvature. 

It should also be noted that there is currently no universal technique for parametric modeling of GTE blade profiles, 
and the selection of a specific technique depends both on the goals of modeling and on the characteristics of the specific 
type of engine being designed. However, the forms of input information can vary significantly. It can be a geometric 
model in the form of a file in the format of one of the CAD systems, or just a set of points from a database of profiles 
of the type [16]. 

Thus, an additional requirement for the parametric modeling technique of the blade profile, due to the purpose of the 
model, is versatility with respect to the format of the source data. 

This research was aimed at the creation of the technique for approximating the profile of the GTE blades, free from 
the above disadvantages and convenient for use in the subsequent analysis of the working conditions of special coatings 
on the surface of the blades. This was demonstrated by the example of solving the problem of gas dynamics for three 
different blade models, whose profile was parameterized through the proposed technique. 

Materials and Methods. The authors used an option of parametric models of the blade profile based on the orthogonal 
Legendre polynomials [17]. A Python application was created for automated processing of source profiles. This made it 
possible to calculate the coefficients of approximating polynomials for the contour lines of the blade, visualize the 
calculation results in the form of graphs of initial points and approximation curves for the suction and pressure sides of 
the blade, save the initial points of the profile, an array of coefficients of approximating polynomials and the coefficient 
of determination of the model in the database, and create a dxf file based on the points of approximating functions to 
transfer it to the CAD system and an engineering analysis package. 

The constructed two-dimensional geometric models of the blades were used to solve the problem of a stationary 
aerodynamic flow around the blade in the approximation of potential flow. The velocity distributions calculated as a result 
of solving this problem in the flow were used to study the effect of the blade profile on its cooling in the aerodynamic flow. 

Calculation method. The orthogonal polynomial systems used as the basis of the proposed method provide linear 
independence of the approximation coefficients and are devoid of a disadvantage that causes difficulties in constructing 
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models that are continuous in parameters. The degree of polynomials is selected to be high enough — of the ninth order, 
so that the set of two polynomials for two parts of the profile (suction and pressure sides) is twenty parameters, which is 
sufficient for approximating very complex profiles. The type of the approximating function and the method of calculating 
estimates of its coefficients are described below. 

We denote the approximating line of the blade profile using function (x). The expression of this line in the case of a 
number of Legendre polynomials of the ninth degree is represented as: 


y(x)=)o ails(x), (1) 


where a; — coefficients that are parameters of the model; L(x) — Legendre polynomials calculated by formulas: 


L, =(3x?-1)/2; 
L =(5x?-3x) /2; 
Ly = (35x* -30x? +3)/8; 
is = (63x3 - 70x? +15x)/8; (2) 
Lo = 231/16x° —315/16x* +105/16x? — 5/16; 
Ly = 429/16x7 —693/16x° +315/16x? —35/16x; 
Lg = 6435/128x* —3003/32x° +3465/64x* —315/32x7+35/128; 


12155 6435 9009 1155, 315 
128x° 32x? 64x° 32x? 128x- 


Let the coordinates of the initial blade profile be given by a set of points yf = y°(x;),i=1,...,N. Substituting 


values x; in functions (2), we obtain matrix X of size N x 10. The coefficients of the approximating polynomial for the 
given blade profiles were found by the least squares method from the minimum condition of the sum of squared deviations 
of the given values yf and calculated through the regression equation (1): 


A N c 2 i 
a=) (9 y(ai, x:)) > min. (3) 
The estimates of coefficients a; were calculated using the well-known regression analysis formula: 
-1 
a=(X"X) XTy, (4) 
where a — notation for the vector of coefficient estimates; y — for the vector of points of given profile y/. 


The approximation accuracy was estimated by the determinacy coefficient of the model (the coefficient of 
determination R?), which, with a good approximation, should be close to unity: 


way Dlv)-yF) 


, (5) 


where y° — average value of y/. 

In design practice, it is customary to describe the leading and trailing edges of the contour of the blades with the radii 
of circles. In this regard, the coordinates of the points of the approximated profiles between the arcs of the circles of the 
leading and trailing edges were considered in this paper. For the purposes of approximation, the contour of the blade in 
question was assumed to be horizontal, so that the blade chord line coincided with the abscissa axis of the Cartesian 
coordinate system. To preserve the orthogonality condition of the models, size normalization is mandatory. In this regard, 
all dimensions were normalized by the length of the blade chord, so that the abscissas of the profile points lay in the range 
of 0.1. This provided the orthogonality of functions (2) with coefficients calculated by formulas (4). 

Research Results. Table 1 shows the coefficients of approximating polynomials for blade profiles belonging to 
different types of GTE to study the heat fatigue of a thermal-protective coating consisting of two metal and ceramic layers: 
profile of the compressor blade C8626 [13], profile of the NASA blade of a highly efficient high-pressure GTE [18], and 
profile of the simulator blade used in [19]. 
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Coefficients of approximating polynomials for blade profiles from different sources 


Table 1 


Coefficients [13] [18] 

suction side | pressure side suction side pressure side suction side pressure side 
do 142.4250 —123.0993 —1275.6471 717.8336 —1043.6707 527.0915 
a —374.3446 324.9875 3328.6126 —1868.9670 2719.5055 —1369.9895 
ar 477.9332 417.7186 —4181.1239 2338.1048 —3406.1466 1707.9198 
a3 —444 4554 392.7204 3794.0526 —2108.4202 3077.0694 —1531.8025 
as 324.6819 292.0489 —2684.7257 1479.0172 —2164.2992 1066.1471 
ds —189.2284 174.3673 1502.4170 —818.4219 1201.7203 —583.7501 
a6 86.8260 —82.4103 —654.4877 351.5869 —518.2516 247.3087 
aq —30.0775 29.5457 211.8930 —111.9355 165.6344 —77.2958 
dg 7.1533 —7,3024 —46.0284 23.8450 —35.3931 16.0493 
ay —0.9024 0.9652 5.1021 2.5861 3.8411 —1.6752 
R2 0.9999 0.9993 0.9896 0.9990 0.9982 0.9994 


Figure | shows the starting points of these profiles and the corresponding approximating lines. 
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Fig. 1. Initial points and approximating lines of blade profiles presented in Table 1, 
sources of blade profiles: a — [13]; b — [18]; c — [19] 
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As mentioned above, the aim of parameterization in this paper is to approximate the profile by a continuous function 
followed by modeling the distribution of external parameters over the surface of the blade. To illustrate the possibility of 
solving this problem, we consider a simple two-dimensional model of a blade in a stationary aerodynamic flow, in the 
approximation of potential flow. In this case, the distribution of the velocities of flow v(x, y) is described by the equation: 


ve , (6) 


where function y is found from the solution to the Laplace equation: 
Ay = 0, (7) 
which in this case corresponds to the incompressibility condition (V v= 0) and vortex-free flow (V Xv= 0). 
The boundary condition for function y will be equal-zero velocity in the normal direction of the blade surface, which 
means that w is constant on this surface. 
We denote the profile line of the blade surface by S and search for a solution to equation (7) in region Q, bounded by 
the outer part with respect to S and the inner part with respect to circle C of a sufficiently large diameter, compared to the 


length of the blade chord, so that the flow is homogeneous on surface C. That is, the conditions are set at the boundary of 
the domain 0OQ=SUC: 


WV] 5 =O. W[6 = Vex, (8) 
where Vv... — uniform flow velocity vector. 


After calculating the velocities, the pressure distribution can be found using the Bernoulli equation: 


2 
v 


2 
where p — gas density; B — Bernoulli's constant, which can be set equal to ps + ae where p. — pressure at an 


infinite distance from the blade. 
The variational formulation of problem (7) looks like this: 


[vvvw=o. Vwe H)(Q). (10) 
Q 


The solution to this problem was carried out by the finite element method. For comparison, the calculation was 
performed for three variants of blade profiles, shown in Table 1, with the same chord length I. 

The obtained results were further used in solving the problem of cooling the blade in an aerodynamic flow. Such a task 
makes sense for thermocyclic testing of coatings on simulator blades. In this case, the problem statement looked as follows: 


0,7 -V-(KVT)+v-VT =0, (11) 
where 7 — temperature; t — time; « — thermal diffusivity coefficient. 
Initial temperature distribution: 

T(t =0, x)=Tp(x) (12) 
was assumed to be as follows: inside the blade, a constant temperature was 400°C, ambient temperature was 0°C. The 
boundary conditions corresponded to the absence of heat flow on the outer contour C: 

ar 
on 


= 0, tle cage (13) 


Cywen>O0 
The distribution of flow velocities in problem (11) was calculated when solving problem (10), while the velocity on 


contour C was assumed to be v=10[//s] (length dimension in units of chord length). The thermal conductivity 


coefficients of the blade and gas were assumed to be 0.1 and 0.01 [! hs s|, respectively. 


The program code for solving problems (10), (11) was written in the input language of the universal finite element 
package FreeFem++ [20]. 
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Figure 2 shows the pressure and temperature distributions calculated as a result of solving problems (11), (12), when 
the blade cools down for 25 seconds in the aerodynamic flow for the three studied blade profiles shown in Figure 1. 
Table 2 shows the calculated values of the minimum relative pressures near the blade surfaces and the maximum 
temperatures of the blades. 


°) 


Fig.2. Calculated distributions of pressure (p, left) and temperature (T, right) when modeling the flow and cooling in the flow of 
the blades of the studied profiles (profile designations — as in Fig. 1). Flow vector direction angle with blade chord direction was 10° 


Table 2 


Minimum relative pressure near the surface of the blade 
and maximum temperature of the blade after 25 s cooling 
in the aerodynamic flow 


Parameters Blade profile 
fs [18] [19] 


Tris (t= = 25), 156.9 p52 | 2 | 2834 4 
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Discussion and Conclusion. Judging by the data presented in Table 1, the coefficients of determination R? for all 
studied blade profiles were close to unity, despite the fact that the shapes of the profiles differed significantly. This 
confirms the good quality of the approximation and is also illustrated by the data in Figure 1. Thus, the proposed method 
makes it possible to construct approximating blade profile lines for various types of GTE with high accuracy. 

The paper also illustrates the possibility of using approximating functions under calculating the distribution of 
pressures and temperatures on the surface of the blade on the basis of solution to the problem of gas dynamics. Further, 
this distribution can be used to optimize the structure of special blade coatings, such as heat-protective coatings. The data 
shown in Figure 2 and in Table 2, confirm that this task can also be successfully solved. The patterns of pressure and 
temperature distributions for different blade profiles have common features, but the detailed picture differs. It can be seen, 
specifically, that the blade of the profile [19], which has a larger cross-sectional area compared to the others, is 
characterized by a higher value of relative pressure in the area of the outlet edge and a lower cooling rate. This is also 
illustrated by the data given in Table 2. Since, from the point of view of the strength of the adhesive bond of the coating 
with the substrate, the most dangerous forces are directed normally away from the surface of the blade, the relative 
minimum pressure near the surface is relevant. Its values are presented in Table 2. With account to the calculation 
obtained, it is possible to reasonably approach the design of a special coating for the surface of a blade of a specific 
profile: adjust the materials and thickness of the coating depending on the degree of loading of individual surface areas. 

Note that these calculation results are presented only as an illustration of the possibilities of the proposed profile 
approximation technique. Therefore, the model of gas dynamics used in the work is quite simplified. In a specific task, it is 
necessary to use a more rigorous model with parameters corresponding to the specified type and mode of the GTE operation. 

Thus, the proposed technique for approximating the blade profile of the GTE, based on the use of orthogonal 
polynomials, is a convenient tool that provides automating the development of a geometric model of the blade and 
comparing different types and profiles of blades, solving the corresponding problems of gas dynamics. At the same time, 
for a given blade profile and working conditions of a GTE, it is possible to obtain a distribution of temperatures and forces 
acting on the blade surface, which is required to predict the durability of special coatings [21]. Based on this calculation, 
it is possible to optimize the coating technology by varying its thickness and/or composition, increasing the coefficient of 
durability in the most dangerous areas. 
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